Demonstration of Turnstiles as a Chaotic Ionization Mechanism in Rydberg Atoms 
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We present an experimental and theoretical study of the chaotic ionization of quasi-one- 
dimensional potassium Rydberg wavepackets via a phase-space turnstile mechanism. Turnstiles 
form a general transport mechanism for numerous chaotic systems, and this study explicitly illumi- 
nates their relevance to atomic ionization. We create time-dependent Rydberg wavepackets, subject 
them to alternating applied electric-field "kicks", and measure the electron survival probability. 
Ionization depends not only on the initial electron energy, but also on the phase-space position of 
the electron with respect to the turnstile — that part of the electron packet inside the turnstile 
ionizes after the applied ionization sequence, while that part outside the turnstile does not. The 
survival data thus encode information on the geometry and location of the turnstile, and are in good 
agreement with theoretical predictions. 

PACS numbers: 32.80.Rm, 05.45. Gg, 05.45.Ac, 45.50.Pk 



Chaotic behavior appears in diverse complex systems, 
over an enormous range of physical scales, including the 
formation of weather patterns, mixing of fluids, firing 
of neurons, and transport in the solar system. Among 
these, photoabsorbtion and ionization in atomic gasses 
have proven to be excellent testbeds for both classical 
and quantum chaos. For example, oscillations in the pho- 
toabsorption spectra of atoms in applied fields have been 
intimately linked to chaotic electron orbits [1], and res- 
onant islands have proven to be barriers to microwave 
ionization [5]. More recently, such resonant islands have 
been used to trap, control, and engineer electronic Ry- 
dberg wave packets [3|. Such experiments highlight the 
utility of highly excited Rydberg electrons as high reso- 
lution probes of chaotic phase spaces. 

Whereas the above examples focused on steady state, 
or nearly steady state, dynamics, the present work fo- 
cuses on revealing the time-dependent mechanism under- 
lying the chaotic ionization of an electron wave packet. 
The experimental protocol is based on the theoretical ob- 
servation [3] (see also Ref. [S]) that the ionization mech- 
anism can be explained in terms of a homoclinic tangle 
and its corresponding turnstile [6l [7j- The turnstile is 
a structure within phase space that promotes the elec- 
tron from a bound to an unbound state, and thus serves 
as the critical step in the ionization process. Turnstiles 
have been theoretically applied to chaotic transport in a 
wide variety of physical systems ^7, . Experimental stud- 
ies, however, are significantly more limited, with notable 
examples to chaotic or turbulent fluids [S], including re- 
cent work on Lagrangian coherent structures IHl [TO] j and 
to optical microcavities [TT]. However, to our knowledge, 
their structure has not previously been experimentally 
measured in atomic ionization. We present experimental 
data on the observation of the phase-space turnstile in a 
system consisting of highly excited quasi-one-dimensional 
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FIG. 1: (Color online.) (a) The experimental timing se- 
quence for the applied electric force F{t). The first pulse is 
the focusing kick. The remaining three pulses are the ioniza- 
tion kicks, (b) The pulse sequence used in the ID simulations 
is a square-wave version of the experimental sequence. 

Rydberg atoms [12] exposed to alternating electric-field 
pulses, or "kicks" . The ionization probability depends 
not only on the electron energy, but, crucially, on the 
phase-space position of the electronic state with respect 
to the turnstile at the moment the ionization kicks are 
applied. We demonstrate that our measurements are a 
sensitive probe of the phase-space position and shape of 
the turnstile. 

Experimental protocol and data: Potassium 
atoms are first excited to a high principal quantum num- 
ber {n « 306 or n « 350) quasi-one-dimensional Ryd- 
berg state using the protocol detailed in Ref. [12]. A 
small electric field "focusing" pulse is applied to the 
atom that delivers an impulsive kick to the electron of 
scaled strength Ap = nAp = 0.08 directed toward the 
nucleus (Fig. [T^). (Tildes denote scaled atomic units: 
f = r/n^, t = t/n^^ E = ri^E.) This creates a non- 
stationary electronic wave packet that undergoes strong 
periodic focusing near the outer classical turning point. 
After some time delay td, a sequence of three alternating 
positive and negative ionization kicks is applied (Fig.jl^). 
These are much stronger than the focusing kick, with 
Ap = 0.25,-0.5,0.25. The fraction of Rydberg atoms 
that survive as a function of delay time is measured 
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FIG. 2: (Color online.) Measured survival probabilities as 
a function of delay time td for the values of peak-to-peak 
duration T indicated. Squares denote experimental data; red 
lines are sinusoidal fits 1131. 



using field ionization. Figure [2] shows results for both 
n « 306 (left column) and n « 350 (right column). The 
delay time tfi is increased from ns to 20 ns, and the peak- 
to-peak duration T of the ionization sequence is adjusted 
between 5 ns and 15ns. The kick durations are all fixed 
at 600 ps. 

The data in Fig. [2] show a clear periodic behavior. 
This is elucidated by fitting to a sinusoid P{td) = Pq + 
Asm(2T:t(i/To + (j)) and extracting the fitting parameters 
Pq, a, To, and (j). The first three parameters are recorded 
in Table |l) the last is plotted in Fig.|3] The following fea- 
tures stand out. (i) The oscillation period Tq is nearly 
independent of T and matches the classical Kepler period 
Tk of the original state; for n — 306, Tk ~ 4.35ns; for 
n = 350, Tk — 6.52ns. (ii) The parameter Pq is the sur- 
vival probability averaged over td- It and the oscillation 
amplitude. A, do not depend strongly on T. (iii) The 
phase shift (p varies considerably with T, as seen in the 
shift of the vertical lines in Fig. [2j 

A classical trajectory Monte Carlo (CTMC) simulation 
is sufficient to reproduce the essential experimental sig- 
natures. We use a one-dimensional hydrogenic model |14j 
with the square-wave forcing shown in Fig. [TJd. The fit- 
ting parameters derived from this model are included in 
Table |I] and Fig. |3] under the label "ID" and are in good 
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TABLE I: Fit parameters Po, A, and To from the experiment 
(exp), ID and 3D simulations, and a ID lobe analysis (ID'). 
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FIG. 3: (Color online.) Phase shift as a function of T. Small 
shifts in T are applied to ID and 3D data to visually separate 
data markers [13) . 



agreement with the experimental data. As shown be- 
low, these parameters reflect the underlying phase-space 
geometry. The amplitude A, however, depends more 
strongly on the details of the electronic initial state, and 
the experimental oscillation amplitude is somewhat less 
than the model predictions. As a consistency check the 
results of more involved three-dimensional CTMC simu- 
lations are also included. 

Connection to turnstile geometry: We now 
demonstrate how the experimental data reveal the pres- 
ence of a phase-space turnstile and how the turnstile ge- 
ometry provides a qualitative and quantitative frame- 
work for understanding trends in the data. Consider 
first the electron dynamics when subjected to periodic 
forcing, in which the single square-wave forcing cycle 
of duration T (Fig. [T]) is replaced by a periodic repe- 
tition (Fig. |4^). A stroboscopic picture of the dynamics 
is defined by the Poincare map {r,pr) i— > (r',p^), which 
takes the radial position and momentum of the electron 
at a given time and returns their values one forcing pe- 
riod T later. This map has a fixed point ai r = oo 
4J to which stable and unstable manifolds are attached 
(Fig. [ijs). (A stable/unstable manifold /W^ consists 
of those phase-space points that converge to the fixed 
point in the positive/negative time direction [7].) The 
stable and unstable manifolds together form a homoclinic 
tangle [7, {a "broken separatrix" ) , which defines the in- 
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FIG. 4: (Color online.) Panel (a) shows the periodic forc- 
ing F{t) used to compute the homoclinic tangle in panel (b). 
The stable/unstable manifolds W^/W" are the thick/thin 
red/blue lines. Kicking parameter values are chosen for figure 
clarity and are not the experimental values used in subsequent 
figures. 



ner gray zone (roughly the "bound" electron states) and 
the outer white zone (the "ionized" states.) The tangle 
also defines regions called lobes, which fall into two cate- 
gories, those that govern electron capture (C„) and those 
that govern electron escape/ionization {E„). Under the 
forcing dynamics. En En+i and C„ i-> Cn+i- The 
critical step for ionization is the E^i i— > Eq transition 
as this promotes electron states from bound to ionized. 
An analogous process Cq i— > C_i governs capture, and 
collectively this mechanism is called a phase-space turn- 
stile [5]. Being only interested in ionization here, we call 
E-i the "turnstile lobe". Its size, shape, and position 
govern the electron survival probability [15 . 

Following the focusing kick, the electronic state is no 
longer stationary. The corresponding classical ensemble 
has a narrow energy distribution p{E) (Fig. [g]) centered 
at £ = -0.5 with AE = 0.053. It forms a partially 
localized, or "focused" , ensemble that hugs and moves 
along the E = —0.5 shell. (Unless otherwise noted, the 
remainder uses the n = 350 case.) The observed oscil- 
lations in survival probability (Fig. [2]) can now be un- 
derstood as the ensemble passing in and out of the E^i 
lobe — survival is higher when the bulk is outside 
and smaller when inside. Figure [5] illustrates this with 
snapshots of the ensemble at a sequence of times: b) Im- 
mediately after the focusing kick td = 0; the ensemble 
roughly corresponds to the energy shell shifted slightly 
downward, c) The first survival maximum; the ensemble 
is almost entirely outside E_ i , resulting in negligible ion- 
ization, d) The first survival minimum; the ensemble has 
reflected once off the nucleus and now lies almost entirely 



within E-i, resulting in large ionization, e) The second 
maximum. Since the ensemble trajectories lie close to 
the original energy shell, the survival oscillation period 
To is approximately the Kepler period T^- 

Next consider the average survival probability Pq; 
1 — Pq can be interpreted as the fraction of time a tra- 
jectory spends inside averaged over the focused en- 
semble. To visualize this, it is easiest to work in the 
canonical coordinates E (electron energy — 1/f) 

and i (time to reach the nucleus) . Consider Fig. [6^. The 
bottom left curve is the negative Kepler period, forming 
the left phase-space boundary. It is physically identi- 
fied with the right vertical boundary line t = 0, which 
represents nuclear impact. Time evolution in these coor- 
dinates consists of uniform motion rightward along hor- 
izontal lines; when a trajectory reaches t = 0, it jumps 
back to t = —Tk- The lobe intersects the energy 
shell E = —0.5 in one large segment. (See Fig. ^p, 
T = 5 ns, for intersection segments.) The lobe has two 
"horns", which stretch downward, intersecting the left 
boundary. From here, they re-emerge as two tendrils 
at the right boundary (taking advantage of the periodic 
boundary conditions), which reach up and again inter- 
sect E = —0.5 in several short segments. In energy-time 
coordinates, the fraction of time a trajectory of energy 
E spends inside E_i equals the total length of the inter- 
section between E_i and the line of constant E, divided 
by Tk- Since AE for the focused ensemble is small, we 
need only consider E = —0.5, with Tk = 27r. Thus, 
1 - Po ~ L/2n, where L = 0.308 is the total length of 
the intersection segments in Fig. [6]d, T — 5 ns. Thereby, 
Pq = 0.692, as recorded in Table |I] column ID', and in 
excellent agreement with the experimental value 0.697. 

Previous work |3] showed that as T increases (with 
kick strength and duration fixed), the lobe in Fig. [6^ 
shifts leftward. This is reflected in Fig. ^p, which shows 
the intersections of E^i with the E — —0.5 line for in- 
creasing values of T. In particular, the primary wide 
segment shifts left with increasing T. However, another 
critical effect is the increase in number and length of the 
shorter segments on the right. As above, these segments 
are the intersections with the horns, which have wrapped 
around. Physically, the wide segment contains trajecto- 
ries that strike the nucleus once during the ionization 
sequence, whereas the shorter segments contain trajec- 
tories that strike the nucleus multiple times during ion- 
ization. As more time elapses between ionization pulses, 
there is more time for trajectories to strike the nucleus, 
and the relative importance of these ionization pathways 
increases. By T = 15ns, 24% of the ionizing trajecto- 
ries have multiple nuclear impacts. The corresponding 
Po values, computed from the interval lengths in Fig.|6]D, 
are recorded in Table |l] column ID' for varying T. These 
agree well with the ID and experimental data. Were ion- 
ization via multiple impacts not included, this agreement 
would be notably worse, especially for T = 15ns. 
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FIG. 5: (Color online.) (a) Experimental survival probability for n = 350, T = 5ns [13]. (b)-(d). Distribution of the classical 
ensemble (dots) relative to -E-i, measured at the times indicated. The thick (green) curve is the E — —0.5 energy shell. 
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FIG. 6; (Color online.) Panel (a) shows the _E_i lobe in 
energy-time coordinates. p{E) illustrates the distribution of 
energies in the focused ensemble. Panel (b) shows the in- 
tersection segments of the E-i lobe with the E — —0.5 line 
for different T values. The dashed line is positioned at the 
negative Kepler period —2tt and the rightmost segments are 
repeated on the left taking advantage of the periodic bound- 
ary conditions. The large dots are the average positions of 
the segments. 



The leftward shift of the intervals in Fig. [6b with in- 
creasing T explains the phase shift in Fig. [3] The large 
dots in Fig. are the average positions (?) of the escap- 
ing points (where the average is taken using the repeated 
segments on the left.) To account for the initial position 
of the focused ensemble, the values of (i) are plotted as 
— {i) — 3tt/2 in Fig. [sj where they track the changes in 
phase seen as T increases. Ionizing trajectories that im- 
pact the nucleus multiple times again play an important 
role, as they shift the average dot leftward in Fig. [6}d. 

Conclusions: Present work identifies the turnstile 
lobe as the critical mechanism for promoting bound elec- 
tronic states to ionization, the turnstile geometry provid- 
ing a convenient framework for explaining the experimen- 
tal results. More broadly, the results demonstrate that 
kicked atomic systems provide a convenient laboratory 



to explore the turnstile mechanism, common to a wide 
variety of physical systems exhibiting chaotic transport. 

We thank Jeff Mestayer for initial assistance; K.B. ac- 
knowledges support of the UCOP Presidential Disserta- 
tion Year Fellowship; K.M. acknowledges the support 
of NSF grant PHY-0748828; and F.B.D. acknowledges 
the support from the Robert A. Welch Foundation under 
grant C-0734 and the NSF under grant PHY 0964819. 



[1 
[2 

[3- 

[4: 

[5 

[6; 

[7 
[8 
[9 
[10 

[11 
[12 
[13 

[14 
[15 



D. Kleppner and J. B. Delos, Found. Phys. 31, 593 
(2001). 

P. M. Koch and K. A. H. van Leeuwen, Phys. Rep. 255, 
289 (1995). 

W. Zhao, J. J. Mestayer, J. C. Lancaster, F. B. Dunning, 
C. O. Reinhold, S. Yoshida, and J. Burgdorfer, Phys. 
Rev. Lett. 97, 253003 (2006). 

K. Burke and K. A. MitcheU, Phys. Rev. A 80, 033416 
(2009). 

K. A. MitcheU, J. P. Handley, B. Tighe, A. Flower, and 
J. B. Delos, Phys. Rev. Lett. 92, 073001 (2004). 
R. S. MacKay, J. D. Meiss, and I. C. Percival, Physica D 
13, 55 (1984). 

S. Wiggins, Chaotic Transport in Dynamical Systems 
(Springer- Verlag, New York, 1992). 

T. H. Solomon, S. Tomas, and J. L. Warner, Phys. Rev. 
Lett. 77, 2682 (1996). 

G. A. Voth, G. Haller, and J. P. GoUub, Phys. Rev. Lett. 
88, 254501 (2002). 

S. C. Shadden, J. O. Dabiri, and J. E. Marsden, Physics 
of Fluids 18, 047105 (2006). 

J.-B. Shim, S.-B. Lee, S. W. Kim, S.-Y. Lee, J. Yang, 
S. Moon, J.-H. Lee, and K. An, Phys. Rev. Lett. 100, 
174102 (2008). 

F. B. Dunning, J. J. Mestayer, C. O. Reinhold, 
S. Yoshida, and J. Burgdorfer, J. Phys. B: At. Mol. Opt. 
Phys. 42, 022001 (2009). 

Due to a triggering problem, the phase shown for the 
n = 350, T = 5ns data was derived from a different 
run, using 67% smaller kick strengths; this is justified by 
the relatively weak dependence of the phase on the kick 
strength. 

Our model ignores core electrons, due to the negligible 
core scattering of the quasi- ID electronic state [12| . 
A slight difference exists between the turnstile defined 
for periodic forcing (Fig. |4]| and that defined for a single 
forcing cycle as in the experiment (Fig. [m . 



